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Abstract. In this paper we present a new, accurate determination of the three components of the absolute space 
velocity of the Galactic globular cluster NGC 6397 (I ~ 338°, b ~ -12°). We used three HST/WFPC2 fields with 
multi-epoch observations to obtain astrometric measurements of objects in three different fields in this cluster. Fhe 
identification of 33 background galaxies with sharp nuclei allowed us to determine an absolute reference point and 
measure the absolute proper motion of the cluster. Fhe third component has been obtained from radial velocities 
measured on spectra from the multi-fiber spectrograph FLAMES at VT2-VLT. We find (fi a cos5, /i<s)j20oo.o = 
(+3.39F0.15, -17.55F0.15) mas yr" 1 , and Kad = +18.36F0.09 (±0.10) km s~\ Assuming a Galactic potential, 
we calculate the cluster orbit for various assumed distances, and briefly discuss the implications. 
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1. Introduction 

After 15 years of activity, the Hubble Space Telescope 
(HST) is producing large volumes of accurate astrometry 
from the measurements that the high spatial resolution of 
WFPC2 and ACS images allows. In the archive are nu- 
merous fields with multi-epoch observations that are well 
separated in time, and first-epoch observations exist for 
even more fields. Fhese are a gold mine for proper-motion 
measurements in a large variety of Galactic and extra- 
galactic fields, including the crowded regions of globular 
clusters (GCs). 

Send offprint requests to: L. R. Bedin 

* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Felescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS 
5-26555, and with the ESO Very Large Telescope + FLAMES 
at the Paranal Observatory, Chile, under the program 71. D- 
0076 



Globular clusters are important Galactic tracers; their 
orbital motion offers insight into the kinematics and dy- 
namics of the halo, and their motion provides important 
clues to the distribution of mass in the Galaxy. Moreover, 
knowledge of GC orbits enables us to study their interac- 
tion with the Galactic disk and halo, and is important to 
dynamical studies of the relationship of the present-day 
mass function of a cluster to the initial one. 

In the first paper of this series (Bedin et al. 2003, Paper 
I) we presented the absolute proper motion of the globu- 
lar cluster M4, using as reference a background QSO. In 
this paper we use the background galaxies in the field of 
another nearby GC, NGC 6397. We make use of archival 
material, as well as of new images taken as a second epoch 
(GO-8656). We also use ground-based spectra from the 
multi-fiber spectrograph FLAMES (Pasquini et al. 2002) 
on VT2-VLT to measure the radial velocity of the clus- 
ter. Finally, with the new measurements of proper motion 
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Table 1. Description of the data set used in this work. 



ID 


date 


EXP- TIME 


FILT 


GO 


first epoch 


A 


Apr'97 


9x ~800s 


F814W 


6797 


B 


Mar'94 


3x ~1000s+2100s+500s 


F814W 


5092 




Mar'94 


5x ~1000s+400s 


F606W 


5092 


C 


Oct'96 


3000s+2x 1900s 


F814W 


6802 




Oct '96 


1300s+700s 


F814W 


7203 




Oct'96 


1700s+400s 


F606W 


6802 




Oct'96 


300s 


F606W 


7203 


D 


Apr'95 


1800s 


F814W 


5369 




Apr'95 


2 x 1200s 


F814W 


5370 


later epochs 


A 


Mar'03 


500s+260s 


F606W 


9676 


B 


Apr'01 


3xl300s+1200s 


F814W 


8656 




Apr'01 


4x 1300s 


F606W 


8656 


C 


Oct'00 


1200s+7x 1300s 


F606W 


8656 




Oct'01 


1200s+5x 1300s 


F814W 


8656 


D 


Apr'01 


2 x 1300s 


F814W 


8656 



Fig. 1. Finding chart of the WFPC2 fields considered for 
this work. The region is 20' x 20' wide, and is centered on 
the center of NGC 6397. The fields measured in this work 
are labeled 'B','C' and 'D'. For completeness we also show 
field 'A', which we did not use here (see text) but will use 
in forthcoming papers. Fields A, B, C, D are centered 5', 
4f5, 8', and 8' from the cluster center, respectively. 

and radial velocity, and an assumed Galactic potential, we 
explore the orbit of the cluster. 

2. Observations, data reduction 

2.1. WFPC2/HST: Proper motions 

In Fig.^we show, in equatorial coordinates, the footprints 
of the three WFPC2 fields of NGC 6397 that we used in 
the present work. The center of the figure coincides with 
the cluster center. In the following we will refer to each 
field as B, C, and D, as labeled in the figure. Field A 
was not used in the present work because the two epochs 
were oriented differently, and this would have required a 
different approach to the reduction, since the chips can no 
longer be treated independently. 

Tabled lists for our three fields the exposures that we 
used. Almost all of the image sets are well dithered, with 
both large and fractional pixel offsets, allowing precise as- 
trometric measurements and a more accurate assessment 
of the errors. 

Positions were measured for each filter, chip, and 
epoch independently, by using the effective-point-spread- 
function (ePSF) algorithms described by Anderson & 
King (2000). As in Paper I, for proper- motion measure- 
ments we used local transformations, and the best distor- 
tion corrections (Anderson & King 2003), after removing 
the 34th-row error (Anderson & King 1999). The photom- 
etry was calibrated to the Vega-mag flight system, follow- 
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Fig. 2. Vector-point diagram of the proper motions in 
equatorial coordinates, combining all the fields. Open cir- 
cles with error bars represent the selected reference galax- 
ies. Dots represent the stars; a circle at 3 mas/yr sepa- 
rates cluster members from field objects. The arrow is the 
adopted absolute proper motion of NGC 6397. 

ing the recipe of Holtzman et al. (1995). For stars with 
>100 DN in their brightest pixel, we were able to measure 
positions in a single image with an error of <0.04 pixel in 
each coordinate. 

In order to avoid introducing systematic errors, we 
measured the proper motions from two epochs in the same 
filter. When two epochs were available for each of two fil- 
ters (fields B and C), the proper motions were obtained 
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from the weighted mean of the two measurements. For 
field B, we linked F814W (2001) with F814W (1994), and 
also F606W (2001) with F606W (1996); in field C we com- 
bined F814W (2001) with F814W (1996), and also F606W 

(2000) with F606W (1996). In field D we matched F814W 

(2001) with F814W (1995). 

A visual inspection of the images reveals many back- 
ground galaxies. Since 33 of these have a point-like nu- 
cleus, we used our ePSF-fitting procedure to measure their 
positions. As was to be expected, the errors found from 
multiple measurements of galaxy positions are much larger 
than those that are typical for stars, and are a strong 
function of the galaxy morphology. Nevertheless, those 
galaxies that give good positions provide an ideal reference 
point for the measurement of absolute proper motions of 
stars. 

Figure shows the vector point diagrams for the three 
fields B, C, and D combined. Small points show the proper 
motions of the stars, and open circles show the motions 
measured for the selected galaxies. For each galaxy we also 
show the proper-motion error bar. The axes are parallel 
to right ascension and declination. Cluster members were 
used to compute the transformations and to derive the 
relative proper motions (see Paper I for details). For the 
pixel scale we adopted the number given by Anderson & 
King (2003), and for the orientations we used information 
in the HST image headers. Cluster members constitute 
the concentrated distribution of points at the lower left. 
The separation between field stars and cluster stars is very 
well defined. Finally, the zero point of the figure was placed 
at the weighted mean motion of the galaxies; this converts 
all the motions to absolute. The centroid proper motions 
of the cluster stars are in good agreement for the three 
fields, as is shown in Table [3 

Taking the weighted mean of the positions of the whole 
sample of 33 galaxies as zero point of the motion, we find 
for NGC 6397 an absolute proper motion, in the J2000.0 
system, of 

fi a cos S = +3.39 ± 0.15 mas yr~ , 

fig = —17.55 ± 0.15 mas yr -1 . 

In Galactic coordinates this corresponds to 

/if cos b = —13.52 ± 0.15 mas yr -1 , 

fi b = -11.70 ± 0.15 mas yr" 1 . 

The dominant source of error is the uncertainty in the zero 
point given by the galaxies. 

We note that the uncertainties of the results for the 
individual fields are quite different. This rests in turn on 
the wide range of the measurement errors for the indi- 
vidual galaxies, as indicated by the lengths of the error 
bars in Fig. 2. Nearly all of the weight comes from the few 
galaxies that have the best measurements; the majority of 
these are in Field B. This shows the importance of finding 
well-exposed galaxies with sharp centers. 

These values are in marginal agreement with those 
found by Cudworth & Hanson (1993), which give a proper 



Table 2. Summary of absolute motion (in [mas/yr]) for 
each field. The number of background galaxies used to 
compute the reference point is also indicated. 



field 


fla cos 5 


C/i a cos S 


1^6 






B 


+3.33 


0.17 


-17.53 


0.17 


12 


C 


+3.14 


0.73 


-17.18 


0.73 


8 


D 


+3.75 


0.40 


-17.81 


0.40 


13 


average 


+3.39 


0.15 


-17.55 


0.15 


33 



motion that is ^13% smaller. However, their values were 
based on a correction from relative to absolute proper 
motion that was obtained from a Galactic model, and 
their procedure seems to provide an analogous discrep- 
ancy (^20% smaller) in the case of M4 (Cudworth & Rees 
1990), when compared with more recent direct determina- 
tions of its absolute proper motion (Dinescu, Girard & van 
Altena, 1999; and Paper I). 

2.2. FLAMES+GIRAFFE/VLT: Radial Velocities 

With the aim of measuring the dispersion of radial ve- 
locities in NGC 6397, in order to get an estimate of the 
distance of the cluster by comparing with the dispersion 
of proper motions (see Bedin et al. 2003b for a brief de- 
scription of the method), we collected spectra of 1508 
stars at the VT2-VLT telescope in May 2003 with the 
FLAMES/GIRAFFE spectrograph (proposal No. 71. D- 
0076). In the present paper we use these spectra to mea- 
sure the mean radial velocity of the cluster. Observations 
were performed in the MEDUSA mode with simultane- 
ous calibration, which allows obtaining about 130 spec- 
tra simultaneously of objects over a field of view of 25 
arcmin in diameter. The high-resolution set-up H9 (now 
called H9B) was used; this set-up allows a spectral resolu- 
tion R = 26, 000 in the 5140-5350 A range. This spectral 
region is well suited for radial-velocity measurements be- 
cause of the large number of metallic lines and the absence 
of H features. It also includes the Mg triplet, which is par- 
ticularly important in radial-velocity determinations for 

metal-poor stars like those of NGC 6397 ([Fc/H] 2.0), 

because it is the strongest feature well visible also in low 
S/N spectra. 

The target stars for radial velocities were selected pri- 
marily on the basis of their position in the color magni- 
tude diagram (CMD). The primary criterion used was to 
select stars along the red giant branch (RGB), the sub- 
giant branch (SGB), and the main sequence (MS) down 
to V = 17.5, on a narrow stripe 0.08 magnitude wide in 
color. This allowed us to limit drastically the contamina- 
tion from field stars. Around the turn-off the stars of NGC 
6397 have a color that is 0.2 magnitude bluer than that 
of the contaminating stars of the Galactic disk. Moreover, 
the evolved stars in NGC 6397 occupy a position in the 
CMD that is well above the location of the nearby disk 
stars. By analyzing in an outer field, far from the cluster 
center, a stripe in the CMD of similar size to our selection 
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strip, we estimated that on average only 12 stars should 
be non-members, i.e., a contamination smaller than 1%. 
(Furthermore, our radial velocities showed that stars in 
this part of the CMD are non- members.) This left virtu- 
ally no contamination in the the sample. 

Data were automatically pre-reduced using the 
GIRAFFE pipeline GIRBLDRS 1 , in which the spectra 
have been de-biased, flat-field corrected, extracted, and 
wavelength-calibrated, using both prior and simultaneous 
calibration-lamp spectra. The resulting spectra have a dis- 
persion of 0.05 A/pixel. A sky correction was applied to 
each stellar spectrum by using the IRAF sarith subroutine 
for subtraction of the average of the sky spectra that were 
observed simultaneously (same FLAMES plate) with the 
stars. 

Radial velocities were obtained from the IRAF- 
fxcor cross-correlation subroutine. Stellar spectra were 
cross-correlated with a synthetic template calculated by 
SPECTRUM 2 , for the mean temperature, log<?, and 
metallicity appropriate for each observed star. The ac- 
curacy of this template was tested on a solar spectrum 
obtained with the UVES spectrograph; the systematic er- 
ror is less than 50 m s _1 . 

From multiple observations of 81 stars it was also pos- 
sible to make a direct estimate of the errors. For stars of 
magnitude V =14, 15, 16, 17, we find errors of 0.2, 0.3, 
0.7, and 1.4 km s _1 , respectively. 

We also searched for possible systematic errors due 
to misalignment of the fibers or to the difference in the 
illumination of the simultaneous calibration fibers and 
the fibers used for the stars 3 . To this end we reduced the 
daytime wavelength-calibration plates (the plates where 
all fibers are pointed at the wavelength-calibration ThAr 
lamp) using the same procedure that we used for our 
stellar plates. Then we measured the velocity shift of 
the fibers, using as a template spectrum the same ThAr 
template that was used for the wavelength calibration. As 
expected, the average radial velocity obtained from the 
spectra from the five simultaneus calibration fibers is zero, 
within the errors. By contrast, the radial velocities from 
the spectra from the remaining fibers show a systematic 
shift, different for each fiber, probably due to their 
incorrect alignment at the entrance of the spectrograph. 
For the purposes of the present paper, it is sufficient 
to consider the mean systematic shift in radial velocity. 
We obtained for the two positioner plates Medusa 1 and 
Medusa 2: 

AV rad [Medusa 1,NGC6397] = -0.175±0.009 km s" 1 , 



1 Blechaet al. (2000). See http://girbldrs.sourceforge.net/ for 
GIRAFFE pipeline, software and documentation. 

2 The LTE spectral synthesis program freely distributed by 
Richard O. Gray. Program and documentation available at 
www.phys.appstate.edu/spectrum/spectrum.html. 

3 See the GIRAFFE pipeline FAQ web page, 
://girbldrs. sourceforge.net/FAQ/index.cgi 



AV rad [Medusa 2,NGC6397] = -0.225 ±0.008 km s" 1 , 

and the mean value is 

AK a d = -0.200 ± 0.025 km s" 1 , 
which must be subtracted from the cluster mean radial 
velocity. 

As a further test we measured the systematic error 
affecting similar observations of another cluster, M4, 
which were collected a month later. In this way we can 
test the time stability of the spectrograph by seeing 
whether these numbers remain constant. We obtained 
in an analogous way for the data-set of M4 the mean 
systematic shift in radial velocity, for the two positioner 
plates Medusa 1 and Medusa 2, 

AFrad [Medusa 1, M4] = -0.176 ± 0.008 km s" 1 

AVrad [Medusa 2, M4] = -0.236 ± 0.008 km s _1 . 

These values agree very well with the previous ones. 
From these results we can assert that GIRAFFE is stable, 
at least for a period of about a month. 

The spatial distribution of our stars, and their helio- 
centric radial velocities as a function of visual magnitude, 
are shown in Figs. [3] and 
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Fig. 3. Spatial distribution of stars with radial velocity 
measurement. The circles indicate the radial distances of 
5 and 10 arc min from the center of the cluster. 

The unweighted average radial velocity 4 , corrected for 
the systematic shift previously discussed, and clipped for 

4 Note that the choice of weights and the calculation of the 
uncertainty in the mean velocity, depend both on the velocity 
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Fig. 4. Radial velocities of the observed sample of stars, 
as a function of magnitude. 

outliers more than 15 km s _1 from the mean, is 

V; ad = 18.36 ± 0.09(±0.03)km s -1 . 

The two errors that are given are first the random error 
and then the uncertainty in the correction for the sys- 
tematic shift. This mean comes from 1486 stars (22 were 
rejected). 

However, there are other effects of the same size that 
afflict spectroscopic measurements of radial velocities, 
such as blue-shift of convective bubbles in giant stars, or 
gravitational redshift in dwarf stars (Dravins, Lindegren 
& Madsen 1999). Since we have observed both giant and 
dwarf stars, we add a global systematic error of 0.10 km 
s^ 1 to our uncertainties. Finally, 

Kad = 18.36 ± 0.09(±0.10)km s -1 . 

[When we analyze the internal dynamics of the cluster, 
we will carefully correct for all these effects, and also each 
fiber will be corrected individually for the calibration bias 
described above.] 

Our heliocentric average radial velocity for NGC 6397 
is within ~ 1.3<7 of the Da Costa et al. (1977) value of 
K-ad = 19.8 ± 1.1 km s _1 , within ~ 2.1<r of the Dubath 
et al. (1997) value of V Ia d = 15.1 ± 1.6 km s _1 , but there 
is large disagreement (~ 5<r) with the K- a d = 18.80 ± 0.10 
km s _1 given by Meylan & Mayor (1991). 

Dubath et al. measured the radial velocity on 
an integrated-light echelle spectrum of the cluster. 

dispersion of the cluster and on the uncertainties of the in- 
dividual velocities (see Pryor & Meylan 1993). A much more 
careful analysis will be done in our future work on the internal 
dynamics. 



Kinematical measurements on this kind of spectra are 
known to be affected by strong observational biases due to 
small-number statistics and strong luminosity weighting 
(of a few red giants). As clearly demonstrated by Zaggia et 
al. (1992, 1993) these biases influence both velocity disper- 
sion and radial velocity: systematic offsets can be as large 
as 1 x the velocity dispersion itself. Therefore it is not sur- 
prising to find a discrepancy between the integrated-light 
radial velocity and the average radial velocity of individual 
stars. 

It is more problematic to explain the discrepancy with 
the Meylan & Mayor (1991) radial velocity, considering 
that errors in their single velocities are similar to ours, 
even though their sample is only ~ 1/12 the size of ours. 
Their measure is based on CORAVEL spectrometer obser- 
vations of 127 RGB and SGB stars. A possible source of 
systematic error could possibly be in a spectral mismatch 
of the CORAVEL mask with the low-metallicity spectra 
of NGC 6397 stars, as Meylan & Mayor indeed mention 
in their paper. Another source of systematic error could 
be the presence of unrecognized field stars, which, consid- 
ering the low velocity of NGC6397 could bias the mean 
radial velocity. Meylan & Mayor observed mainly RGB 
stars in NGC 6397, while the majority of our targets are 
turn-off and MS stars, which should be relatively free of 
field contamination. Unfortunately a direct comparison of 
the radial velocities of the two datasets is not possible. 

In conclusion, considering the work done to limit all 
possible systematics in the GIRAFFE spectra, we believe 
that with a final systematic error < 0.1 km s -1 and the 
large-number statistics of nearly 1600 radial velocities, our 
mean radial velocity is a significant improvement over pre- 
vious measurements. 



Table 3. Input conditions for orbit calculation. Space ve- 
locities in km s^ 1 . 





Distance 1.8 kpc 






u 


V W = Z 


n 


e 


34 ±5 


-92 ±11 -94 ±10 


22 ±5 


131 ± 11 




Distance 2.0 kpc 






U 


V w = z 


n 


e 


41 ±5 


-103 ±11 -105 ±10 


28 ±5 


121 ± 11 




Distance 2.2 kpc 






U 


V w = z 


n 


e 


48 ±5 


-114 ±11 -116±10 


35 ±5 


111 ± 11 




Distance 2.4 kpc 






U 


V W = Z 


n 


e 


55 ±5 


-125 ±11 -126 ±10 


42 ±5 


101 ± 11 



3. Calculation of the orbital parameters 

Our cluster radial velocity and absolute proper motion 
allow us to derive the three velocity components of NGC 
6397, and estimate its Galactic orbit. The new orbit will 
help us, in principle, to trace back the dynamical history 
of the cluster, clarify how close it goes to the Galactic 
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X [kpc] R [kpc] 

M 1 1 i 4 n 



Table 4. Orbit parameters, for the four different dis- 
tances. Units: rf[kpc], L z [kpc km s _1 ], .Etot[10 x km 2 s -2 ], 
P[Myr], i? a [kpc], i? p [kpc], z max [kpc], e[pure number]. 



d 




-Etot 


P 


Ra 


Bp 


^max 


e 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1.8 


933 


-10203 


112 


6.92 


3.34 


1.84 


0.35 


2.0 


846 


-9899 


102 


6.78 


2.94 


2.02 


0.39 


2.2 


766 


-9634 


99 


6.74 


2.83 


2.80 


0.41 


2.4 


693 


-9424 


93 


6.59 


2.57 


3.11 


0.44 





d = 2200 pc 



R [kpc 









: d = 


200C 


pc : 

i , , , ■ 




Fig. 5. Calculated orbits, for assumed distances 1.8, 2.0, 
2.2, and 2.4 kpc. 



center, and assess the possible impact of the motion on its 
internal dynamics and its mass function. 

Orbit integration for a globular cluster requires adopt- 
ing a model of the potential of the Milky Way. We chose 
that of Allen & Santillan (1991), which was constructed 
to fit an assumed rotation curve and a distance and ve- 
locity of the Sun with respect to the Galactic Center. The 
potential is time-independent — clearly a crude approxi- 
mation, because a significant variation of the Galactic po- 
tential is expected over a typical globular-cluster lifetime. 
Nevertheless, it is reasonable to believe that the Galactic 



potential has not changed much in the last few Gyr, so 
that the derived parameters for the orbit of NGC 6397, 
such as the apo- and perigalacticon, can be considered to 
be reasonable estimates. 

The Galactic potential of Allen & Santillan (1991) as- 
sumes density distributions in bulge, disk, and halo com- 
ponents. This gravitational potential is time- independent, 
axisymmetric, fully analytic, and mathematically very 
simple. It has already been used to derive the Galactic 
orbits of nearby stars (Allen & Martos 1986), open clus- 
ters (Carraro & Chiosi 1994) and disk and halo GCs 
(Odcnkirchen & Brosche 1992). The initial conditions are 
given in Table [3] for different cluster distances, and as- 
suming R = 8 kpc, and the value of V /R derived in 
Paper I. 

The integration routine is a fifteenth-order symmet- 
ric, simplectic Runge-Kutta, using the Radau scheme 
(Everhart 1985); it guarantees energy and momentum con- 
servation at the level of 10 -12 and 10 -9 , respectively, over 
the whole orbit integration. The orbits have been inte- 
grated back in time for 5 Gyr, and are shown in Fig. 

The orbital parameters are summarized in Table 0] 
where column (1) lists the adopted cluster heliocentric 
distance, column (2) the z-component of the angular mo- 
mentum, column (3) the total energy, column (4) the or- 
bital period (for radial motion and for vertical motion), 
columns (5) and (6) the apo- and peri-center of the or- 
bit, column (7) the maximum vertical distance the clus- 
ter reaches and column (8) the eccentricity, defined as 
(R a — Rp)/{Ra + Rp)- The orbit parameters that we found 
do not differ much from those published by Dauphole et al. 
(1996) and by Dinescu et al. (1999), which means basically 
that the Galactic model we adopt in this paper is con- 
sistent with those used by the previous authors. Indeed, 
Dauphole et al. (1996) used an earlier version (Allen & 
Martos 1986) of the Galactic model we use here (Allen 
& Santillan 1991). In any case, our newly derived space 
velocities make our results an improvement over previ- 
ous investigations. In particular, this is the first time that 
proper motions calibrated directly against extragalactic 
sources have been used in studying the orbit of NGC 6397. 

Independent of the initial conditions chosen, the orbit 
of NGC 6397 has a boxy nature. The cluster oscillates 
rapidly through a dense part of the Galactic plane, only a 
few kiloparsecs from the Galactic center. This kind of orbit 
makes NGC 6397 very vulnerable to tidal shocks, with a 
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consequent preferential loss of low mass stars, which, as 
first recognized by Piotto, Cool, & King (1997, see their 
detailed discussion), is the likely cause of its anomalously 
flat mass function. 

As the heliocentric distance increases (from bottom to 
top panel) the cluster tends to have a longer period, to 
reach greater heights above the Galactic plane, and to 
show a larger epicyclic amplitude, dipping closer to the 
Galactic Center. The assumed distance has a complicated 
effect on tidal shocks, which work over time to flatten 
the mass function of the cluster. Tidal shocks occur three 
times in each orbital period: at perigalactic passage and 
at each passage through the Galactic plane. For a greater 
assumed distance the cluster is shocked more frequently, 
and the shocks due to the Galactic center are stronger, 
but the shocks due to the plane occur at higher cluster 
velocity and are thus less effective. 

A more accurate distance determination is in progress, 
based on the same observing material that is presented in 
this paper, ft will be a fundamental step in our under- 
standing of the motion of NGC 6397. 
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